Theoretical studies of equilibrium beta limit in heliotron plasmas

In this study, the MHD equilibrium limit of LHD plasmas is studied using a 3D MHD equilibrium calculation code, HINT2, which is an initial value solver based on the relaxation method without the assumption of nested flux surfaces. For finite beta equilibria, flux surfaces become stochastic in the peripheral region. However, though the axis shifts until the conventional limit, the separatrix to limit the equilibrium does not appear. For high beta equilibria, the force balance starts breaking in the stochastic region. To keep the force balance, the pressure gradient in the stochastic region decreases and the fixed profile is reduced. As the result, the volume averaged beta saturates. The beta value, where the force balance starts breaking inside the stochastic region, is proposed as an index of equilibrium beta limit in heliotron plasmas.